function output = fn_rho1(s_grid,phi)

% output = phi(1) + phi(2).*(s_grid-1);

output = 1 + phi.*(s_grid-1);

end